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ABSTRACT 

We present new measurements of the large-scale bulk flows of galaxy clusters 
based on 5-year WMAP data and a significantly expanded X-ray cluster cata- 
logue. Our method probes the flow via measurements of the kinematic Sunyaev- 
Zeldovich (SZ) effect produced by the hot gas in moving clusters. It computes 
the dipole in the cosmic microwave background (CMB) data at cluster pixels, 
which preserves the SZ component while integrating down other contributions. 
Our improved catalog of over 1,000 clusters enables us to further investigate pos- 
sible systematic effects and, thanks to a higher median cluster redshift, allows us 
to measure the bulk flow to larger scales. We present a corrected error treatment 
and demonstrate that the more X-ray luminous clusters, while fewer in number, 
have much larger optical depth, resulting in a higher dipole and thus a more 
accurate flow measurement. This results in the observed correlation of the dipole 
derived at the aperture of zero monopole with the monopole measured over the 
cluster central regions. This correlation is expected if the dipole is produced by 
the SZ effect and cannot be caused by unidentified systematics (or primary cos- 
mic microwave background anisotropies) . We measure that the flow is consistent 
with approximately constant velocity out to at least ~800 Mpc. The significance 
of the measured signal peaks around 500 h^ c fMpc, most likely because the con- 
tribution from more distant clusters becomes progressively more diluted by the 
WMAP beam. We can, however, at present not rule out either that these more 
distant clusters simply contribute less to the overall motion. 
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Subject headings: Cosmology - cosmic microwave background - observations - 
diffuse radiation - early Universe 

The large-scale isotropy of the Universe and the small-scale inhomogeneities that evolved 
into galaxies are thought to originate during inflationary expansion in the early Universe. 
Inflation posits that the primeval space-time was inhomogeneous and this structure should 
have been preserved on sufficiently large scales. At minimum large-scale peculiar velocities 
should arise from gravitational instability caused by mass inhomogeneities seeded during 
the inflationary expansion. On scales > lOOMpc the standard inflationary scenario leads to 
robust predictions for these velocities. There the initial Harrison-Zeldovich slope of the mass 
fluctuations is preserved, so peculiar velocities induced by gravitational instability must 
decrease linearly with scale (e.g. Kashlinsky & Jones 1991); for the concordance ACDM 
model V ims ~ 250( looft ~ lMpc ) km/sec at rf>50-100/i~ 1 Mpc. 

Our discovery of a coherent large-scale flow of galaxy clusters with significantly larger 
amplitude than expected out to ~ 300Mpc (Kashlinsky et al 2008, 2009 - KABKE1,2) 
represents a challenge to the gravitational instability paradigm. Such a "dark flow" could 
indicate a tilt created by the pre-inflationary inhomogeneous structure of space-time (Turner 
1991, Grishchuk 1992, Kashlinsky et al 1994, KABKE1) and might provide an indirect probe 
of the Multiverse. Various explanations have been put forward, including that the flow 
points to a higher-dimensional structure of gravity (Afshordi et al 2009, Khoury & Wyman 
2009), or that it reflects the pre-inflationary landscape produced by certain variants of string 
cosmology (Mersini-Houghton & Holman 2009, Carrol et al 2008). 

Making use of an expanded cluster catalog and deeper WMAP observations, we have 
worked to verify, and expand the Dark Flow study through a program we have dubbed 
SCOUT (Sunyaev-Zel'dovich Cluster Observations as probes of the Universe's Tilt). First 
results from this experiment are reported here. 

1. Data and analysis 

KABKE1,2 and this work utilize a method of Kashlinsky & Atrio-Barandela (2000, here- 
after KA-B), which measures CMB dipole at the locations of X-ray clusters. When averaged 
over many isotropically distributed clusters moving at a significant bulk flow with respect to 
the CMB, the kinematic term dominates the SZ signal, thereby enabling a measurement of 
Ubuik over that distance. In Atrio-Barandela et al (2008, hereafter AKKE) we demonstrated 
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that 1) the thermal SZ (TSZ) signal from clusters extends well beyond the measured X-ray 
extent 0x, and 2) the intra-cluster gas distribution is well approximated by the NFW profile 
(Navarro et al 1996) expected for dark matter in a ACDM model. The temperature Tx of 
hot gas distributed according to these profiles decreases significantly from the cluster cores to 
the cluster outskirts (Komatsu & Seljak 2001), consistent with current measurements (Pratt 
et al 2007) and numerical simulations (e.g. Borgani et al 2004). Consequently, the monopole 
produced by the TSZ component (oc tT x ) decreases, as we increase the cluster aperture, 
whereas the dipole due to the kinematic SZ (KSZ) component (oc r, the optical depth due 
to Thomson scattering) remains measurable out to the aperture where we still detect the 
TSZ decrement in unfiltered maps (KABKE2). As in KABKE1,2 our dipole coefficients 
are normalized such that the dipole power C\ due to a coherent motion at velocity Vbuik is 
Ci, ki n = T 2 MB (r> X 2 ulk /c 2 , where T CMB = 2.725K. We adopt fi total = 1, fi A = 0.7, H = 70h 70 
km/sec/Mpc. 

To improve upon the all-sky cluster catalogue of Kocevski & Ebeling (2006) used by 
KABKE1,2, we have screened the ROSAT Bright-Source Catalogue (Voges et al. 1999) using 
the same X-ray selection criteria (including a nominal flux limit of 1 x 10 -12 erg/sec, 0.1-2.4 
keV) as employed during the Massive Cluster Survey (MACS, Ebeling et al. 2001), as well 
as the same optical follow-up strategy. Unlike MACS, we apply neither a declination nor 
a redshift limit though, thereby creating an all-sky list of cluster candidates that extends 
to three times fainter X-ray fluxes than in KABKE1,2. Optical follow-up observations of 
clusters identified in this manner and lacking spectroscopic redshifts in the literature (and 
MACS) are well underway using telescopes on Mauna-Kea/Hawai'i and La-Silla/Chile. As 
a result, our interim all-sky cluster catalogue currently comprises in excess of 1,400 X-ray 
selected clusters, all of them with spectroscopic redshifts. X-ray properties of all clusters 
(most importantly total luminosities and central electron densities) are computed as before 
(KABKE2). Within the same z-range [z < 0.25) as our previous study, our new catalog 
comprises 1,174 clusters outside the KP0 CMB mask. To eliminate low- mass galaxy groups 
we require that clusters feature Lx > 2 x 10 43 erg/sec; 985 systems meet this criterion. This 
sample represents a significant improvement over the one in KABKE1,2 largely because of 
the substantial increase (z me dian < 0.1 in KABKE1,2 vs z me dian — 0.2 below) in median 
cluster redshift and the higher fraction of intrinsically very X-ray luminous systems (Fig. 
la). 

We applied this new cluster sample to the 5-year WMAP CMB data, processed as 
described in detail in KABKE1,2. The all-sky dipole in the foreground-cleaned maps from 
http: //www.lambda.gs fcTnasa.gov| was removed, and standard CMB masking was applied. 
We also need to remove the primary CMB fluctuations, produced at last scattering, as they 
are highly correlated and would contribute significantly to the measured dipole. To this 



-4- 



end, all maps were filtered as in KABKE1,2 with a filter that minimizes ((ST — 5Tacdm) 2 )- 
The error budget associated with our filtering is discussed by Atrio-Barandela et al (2010, 
hereafter AKEKE). 

Measurement errors were computed as in KABKE1,2 with one important correction. 
Although the filtering removes much of the CMB fluctuations, a residual component remains, 
due to cosmic variance and imperfections of the theoretical model. Since this residual is 
common to all WMAP bands, a component of the errors is correlated between the various 
DA maps. We address this issue in the following manner. For each of the eight differential 
assemblies (DAs) we simulate 4,000 realizations with N c \ =100-1,000 randomly selected 
pseudo-clusters outside of the cluster pixels and the CMB mask. In each realization we 
select the same pseudo-clusters for all DAs and evaluate the mean monopole and dipole 
averaged over all DAs: a , a lm . For the 4,000 realizations at each N c \ we compute the mean 
and dispersion of a , a\ m over all the realizations. The distributions have zero mean and their 

— 1/2 

dispersion gives errors which scale as N cl . We find to good accuracy that the distribution 
of the simulated dipoles (and monopoles) is Gaussian and the errors on each of the averaged 
dipole components are a lm ~ 15A/3/iV c i/iK and on the monopole a ~ 15^/l/N c \jj,K as 
explained in great detail in AKEKE; the errors of the x/z dipole component are slightly 
larger /smaller because of the Galactic mask (Fig. lb). 

2. Results 

We used the filtered maps to compute the dipole and monopole terms, a\ m , a l for each 
DA of the WMAP Q, V, W bands (i = 1, 8) for clusters in cumulative redshift bins up to 
a given z. The results were averaged to obtain the mean values over all eight DAs, ai m , a . 
Since the volume probed out to low z is too small for meaningful measurements, Table 1 lists 
results only for z-bins with sufficient signal-to-noise ratios (S/N). In KABKE1,2 we computed 
the dipole in progressively increasing apertures but no larger than 30' to prevent geometric 
biases from very nearby Coma-type clusters. To further reduce any selection effect related 
to the apparent X-ray extent measured by ROSAT, we here impose a constant aperture for 
all clusters (see Table 1) and compute the dipole component for each z-bin at the constant 
aperture at which the monopole (initially negative because of the TSZ component) vanishes. 

The improved cluster catalog allows us to extend our study to higher z, and to further 
test the impact of systematics. We create Lx-limited subsamples which achieves two impor- 
tant objectives. 1) As the Lx threshold is raised, fewer clusters remain and the statistical 
uncertainty of the dipole increases (oc 1/ \/ N c \ (L > Lx ) ) • If, however, all clusters are part 
of a bulk flow of a given velocity \4>uik, very X-ray luminous clusters will produce a larger 



-5 - 



CMB dipole (oc rVbuik), an effect that might overcome the reduced number statistics, giving 
a higher S/N in the measured dipole. 2) Since, as outlined under (f), the dipole signal 
should increase with cluster luminosity, whereas systematic effects can be expected to be 
independent of Lx, an actual observation of such a correlation would lend strong support to 
the validity of our measurement and the reality of the "dark flow". 

Fig. la shows that the depth to which we probe the flow increases dramatically as 
the Lx-threshold is raised. Table 1 shows the results for each subsample. The monopole is 
strongly negative in the smallest apertures due to the dominance of the TSZ component in 
the central regions and the dipole is shown at the aperture where monopole vanishes. Of the 
three dipole components, the ^-component is best determined, its value always remaining 
negative and its S/N increasing strongly with increasing L X - As shown in Fig. lc the ampli- 
tudes of the latter and of the monopole in the central parts are strongly, and approximately 
linearly, correlated. This correlation provides strong evidence against unknown systematics 
systematics - or primary CMB - causing our measurement. 

Table 1 quantifies our finding of a statistically significant dipole out to the largest scales 
probed (~ SOO/i^Mpc). In KABKE2 we discussed in detail why this dipole is unlikely to be 
produced by systematics; we briefly revisit the issue here: 1) At high statistical significance 
the dipole originates only at cluster positions, and must thus originate from CMB photons 
that passed through the hot intra-cluster gas. For the same reason the dipole can not be 
due to a residual contribution from the all-sky CMB dipole. HEALPix ANAFAST routines 
(Gorski et al 2005), employed in the analysis, further remove any all-sky dipole before the 
filtered maps are produced. 2) Since the dipole is measured at zero monopole, the contribu- 
tions from TSZ and other cluster emissions to the dipole are negligible. 3) The variations 
in the final aperture where the dipole is measured were very small in KABKE1,2, and the 
dipole signal remains in this work which uses a fixed aperture for all clusters. So the dipole 
is not affected by the variations in cluster Q x , which in any case is much smaller than the 
final apertures in KABKE1,2. 4) The measured CMB quadrupole is significantly different 
from that of the ACDM model, so a significant part of the CMB quadrupole is not removed 
by our filter and could leak into other multipoles via the mask. However, we set the filter to 
zero at £ < 4 (KABKE2, AKEKE) for the final maps. More importantly, we have modified 
the pipeline to remove the all-sky quadrupole from the original maps and find no notice- 
able difference in the dipole computed at the cluster locations. This also removes the fully 
relativistic components from the local motion t>i OC ai, down to (fi OC ai/c) 3 corrections to the 
octupole. (5) Finally, we demonstrate that more luminous clusters make a larger, and sta- 
tistically more significant, contribution to the dipole, as is expected if all clusters participate 
in the same flow, independent of Lx- Note also that the intra-cluster medium motions from 
cluster mergers have random directions and thus average down in large catalogs, contributing 
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negligibly to the noise budget below. 

Although the final dipole is measured at zero monopole, tests of cross-talk were con- 
ducted as in KABKE2 by constructing CMB maps from TSZ and KSZ components of vary- 
ing Vb U ik using the derived catalog parameters, but randomly placed clusters. The param- 
eters from randomly placed clusters were compared with those from the original clusters 
- the cross-talk effects are small with results similar to Fig. 6 of KABKE2 (AKEKE). 
To test the robustness further we select 4,000 random subsets of clusters within a given 
configuration and compute the mean dipole and its dispersion. The distribution of the 
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dipoles is Gaussian, the dispersion scales as iV cl and the results are consistent with all 
subsets of clusters moving in the same way within the estimated errors. E.g. randomly 
selecting 250/150 L x > 2 x 10 44 erg/sec clusters out of 322/208 at z < 0.25/0.2 gives 
(a lx , a ly , a lz ) = (3.7 ± 1.8, -4.1 ± 1.7, 4.2 ± 1.5)/(3.6 ±2.1, -5.7 ± 2.3, 4.5 ± 1.9)/iK. 

To calibrate our dipole measurement in terms of an equivalent bulk velocity, we proceed 
as in KABKE2. This still suffers from a systematic bias which overestimates the amplitude 
of the velocity. More importantly, we measure the dipole from the filtered maps, and the 
convolution of the intrinsic KSZ signal with the filter can change the sign of the former for 
NFW clusters. The TSZ signal, being more concentrated as shown in Fig. 9 of KABKE2, is 
less susceptible to this effect. We therefore currently constrain only the axis of the motion; 
the direction along this axis should result from future applications of the KA-B method 
particularly to the 217 GHz Planck data, where the TSZ component vanishes and the angular 
resolution is 5', a good match to the inner parts of clusters at z~0. 1-0.2. Our present pipeline 
computes the cluster properties (central electron density, n e ,o, Tx, core radius -R CO rc) assuming 
a /3-model (/3=2/3). This model has been shown by us to be deficient at the cluster outskirts 
and must be replaced by the NFW profile (AKKE). We hope to accomplish this difficult task 
in the future with better CMB (Planck) and X-ray (Chandra/XMM) data. To summarize, 
our current calibration may overestimate the amplitude of the flow and, strictly speaking, 
we currently measure only the axis of motion. We stress, however, that the existence of the 
flow itself is not affected by this systematic uncertainty. We generated CMB temperatures 
from the KSZ effect for each cluster and estimate the dipole, C^ioo, contributed by each 100 
km/sec of bulk- flow in each Lx, 2-bin. Since the /3-model still gives a fair approximation to 
cluster properties around Qx, we present in Table 1 the final calibration coefficients evaluated 
at apertures of 5' and Qx in radius. When averaged over clusters of all X-ray luminosities 
the mean calibration is a/ (Ci,ioo) — 0.3/iK in each of the z-bins. Within the uncertainties, 
the dependence of the calibration on L x is in good agreement with the measured dipoles, 
particularly for the most accurately measured y-component. Table 1 also shows the mean 
central optical depth, (tq), evaluated from the cluster catalog. Its variation with Lx is also 
in good agreement with that of the measured dipole, which indicates that the clusters can 
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indeed be assumed to have similar profiles. 

With the calibration factors in Table 1, our results are compatible with a consis- 
tently coherent flow at all z. Since the different cluster subsamples probe different depths 
(-^mean/median), we proceed as follows to isolate the overall flow across all available scales: for 
the flow which extends from the smallest to the largest z in each z-bin, we model the dipoles 
as a" m = a n V m , where a n is the calibration constant for clusters in n-th luminosity bin. (Note 
that for each z-bm the luminosity bins are statistically independent). Both a/Ci,ioo an d ( r o) 
scale approximately linearly with the better measured dipole coefficient, ai tV , as they should 
in case of a coherent motion; the linear correlation coefficients are r =0.92/0.93 for corre- 
lation of —d\,y with columns 7/8. We then compute by regression the velocity components 
with their uncertainties using the L^-divisions at each z-bm. These numbers are shown in 
the summary row following each z-bm. quantities for a = a/Cl^loo at 5'. We also computed 
them for a n given by a/C^ioo & t @v and by (r ) in column (7) normalized to the observed 
mean value of (Ci,ioo)- Both give results essentially identical to the ones shown in the Table. 
The latter approximation is equivalent to assuming that all clusters have universal profiles, 
so that the final effective optical depth is oc (r ) x (reduction factor). The results here are 
consistent with KABKEf,2 measurements on smaller scales (^400Mpc), which with revised 
errors become (a lx ,a ly ,a lz ) = (0.7 ± 1.2, -3.3 ± 1.1, 0.5 ± l.)/(0.6± 1.2, -2.7 ± 1.1, 0.6 ± 1.) 
for z < 0.2/0.3 with ^/C^ ~ 0.3/xK (Table 2, KABKE2). 

The dipoles are larger for greater Lx clusters consistent with the dipole originating 
from the bulk motion of the clusters. We note the apparent trend in the central values of 
the better determined y-component peaking at z < 0.16 and decreasing towards higher z. It 
is likely that this decrease is due to dilution of progressively more distant clusters, as shown 
by their smaller monopoles. Nevertheless, it is also possible, in principle, that the flow is 
dominated by the z < 0.16 clusters with the more distant clusters contributing little to the 
dipole. 

3. Discussion 

We find a high likelihood of the existence of a coherent bulk flow extending to at 
least z ~ 0.2 with an amplitude and in a direction which are in good agreement with our 
earlier measurements. Our result constitutes a significant improvement in that it extends 
our previous work to approximately twice the distance accessible to KABKE1,2, supporting 
their hypothesis that the flow likely extends across much (or all) of the Hubble volume. The 
flow's axis is also consistent with earlier measurements of the local cluster dipole (Kocevski 
et al. 2004) as well as with independent measurements of bulk flows on smaller scales by 
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Watkins et al (2009). The velocity reported there is smaller than the numbers in Table 1, 
although the two amplitudes agree at < 2-a level. Agreement between the two sets of central 
values would require a/ (C*i,ioo) ~0.4-0.5/iK, or a reduction by a factor of ~ 2 from unfiltered 
values for NFW cluster profiles. Feldman et al (2009) extend the Watkins et al analysis and 
find that the absence of shear in their flow at <50-100Mpc is consistent with the KABKE1 
suggestion of the attractor at superhorizon distances. 

Fig. [2] displays the results obtained in this study compared to expectations from the 
concordance ACDM model for 95% of cosmic observers. These results cast doubt on the 
notion that gravitational instability from the observed mass distribution is the sole - or even 
dominant - cause of the detected motion. If the current picture is confirmed, it will have 
profound implications for our understanding of the global structure of space-time and our 
Universe's place in it. 

We acknowledge NASA NNG04G089G/09-ADP09-0050 and FIS2006-05319/GR-234 grants 
from Spanish Ministerio de Education y Ciencia/ Junta de Castilla y Leon. 
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(°) d ~ 385 - 755h TO 1 Mpc; (V x , V y ,V z ) = (313 ± 308, -707 ± 265, 643 ± 259) x 
( 6 ) d ~ 385 - 755/if 1 Mpc; (V x ,V y , V z ) = (352 ± 304, -713 ± 262, 652 ± 256) x 



+o.3 M y 
\J (^1,100) 

+0.3M-ff 



km/sec; Veulk = (1, 005 ± 267) X 
km/scc; V Bu lk = (1, 028 ± 265) X 



V (Ci,ioo) 



km/sec; (l ,b ) = (296 ± 29, 39 ± 15)° 
km/sec; (Z , b ) = (296 ± 28, 39 ± 14)° 



Note. - 
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Figure and table captions: 

Table 1: : Columns are: (1) - the limit of the cumulative z-bin. (2) luminosity 
range of the differential Lx-bin. (3) the number of clusters per bin. (4) mean/median 
redshift in the bin. (5) Dipole coefficients, averaged over the eight WMAP DA's over the 
clusters in the bin with 1-a errors, 0\ m . (6) Dipole amplitude, \fC[ = ^ m a\ m . The 
error on the dipole amplitude is derived as a\ = ^ m {d\/C[/ dai m ) 2 a\ m = ^2 m ai m o'l m /C 1 . 
The binning in luminosity was designed, wherever possible, to bin different z-bins by the 
same luminosity range. Note that a ly are always negative. When we select 598 clusters 
out to z < 0.25 with Lx > 10 44 erg/sec the dipoles are consistent with the brightest 
Lx-bin: (di x , di y , di z ) = (—0.3 ± 1.2,-2.5 ± 1.1, 2.4 ± 1.1)//K. (7) Central optical depths, 
tq = y/TT<7Tn e flR COTe , averaged over clusters in the bin derived from our cluster catalog as 
described in the text. (8) Calibration factors, a/Ci,ioo ; f° r clusters in the given bin at 5' 
radial distance from the cluster centers and at 1 X-ray extent for (5- model with /3 — 2/3. (9) 
Measured monopole, over the fixed aperture with the radius shown, after averaging over all 
DA's. 

The row at the end of each z-bin sums up the bulk flow parameters (scale, components, am- 
plitude and direction to the axis of motion) assuming a coherent motion for all the Lx-bins. 
The depth is defined as d = ^mcdianC-f^ 1 . 

W There are only five clusters in this Lx-range at z > 0.12, so for brevity this configuration's 
parameters are not repeated the remaining z-bins, although for completeness it is included 
in bulk flow evaluations ( b \ 

O) ^he ] 3U jj c fl ow j s derived using column 8 parameters without the lowest L^-clusters which 
do not extend beyond z = 0.12 M (see text). 

^ The bulk flow is derived using column 7 parameters including the lowest Lx-clusters ^ 
(see text). 

Figure 1: (a) ^-distribution of clusters in the Lx-bins used in the analysis. Black/blue/green/red 
colors correspond to L x = (0.2 — 0.5, 0.5 — 1, 1 — 2, > 2) x 10 44 erg/sec respectively, (b) Stan- 
dard deviations from simulations of random pseudo-clusters vs N c \. Pluses/ diamonds/triangles 
correspond to a,\ x ja\yja\ z . (c) Correlation between the ^-component of the dipole and the 
central monopole from Table 1 (filled circles). Both quantities increase in amplitude with Lx 
because the optical depth and central temperature increase for more massive clusters: KSZ 
dipole scales as r and the monopole as rTx- The linear correlation coefficient for the circles is 
shown in the upper left. Triangles correspond to differential ^-configuration, 0.12 < z < 0.25, 
where we still have enough luminous clusters for reasonable S/N. Open/filled triangles corre- 
spond to L x > (1,2) x 10 44 erg/sec with 418/260 clusters. When the triangles are included 
the linear correlation coefficient becomes 0.85. 
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Figure 2: Bulk velocity vs depth from Table 1: blue/cyan/grean/red corresponds to 
z < 0.12/0.16/0.2/0.25; parameters for fits (a) are chosen for brevity. Solid/dashed lines 
correspond to the rms bulk velocity for the concordance ACDM model for top-hat/Gaussian 
windows. Black-shaded regions shows the 95% confidence level of the model (see KABKE1 
for details). 



